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Coupling efficiencies of an electromagnetic field with a Kolmogorov phase statistics into a step- 
index fiber in its monomode regime of wavelengths are computed from the overlap integral between 
the phase screens and the far-field of the monomode at infrared wavelengths. 

The phase screens are composed from Karhunen-Loeve basis functions, optionally cutting off 
some of the eigenmodes of largest eigenvalue as if Adaptive Optics had corrected for some of the 
perturbations. 

The examples are given for telescope diameters of 1 and 1.8 m, and Fried parameters of 10 and 
20 cm. The wavelength of the stellar light is in the J, H, or K band of atmospheric transmission, 
where the fiber core diameter is tailored to move the cutoff wavelength of the monomode regime to 
the edges of these bands. 
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I. SCOPE 

One of the promises of monomode fiber optics is the 
removal of corrugations across pupil planes ("cleaning" 
of the beams) by symmetric weighting of the electric 
field across the pupil [3, [l?], [H, |42j]. The transforma- 
tion of phase screens is studied by computation of the 
overlap with the anticipated far-field distribution of the 
fiber optics of the detector system [33[ • We demonstrate 
the negative impact of the growing number of speckles 
[28l | on the throughput of the spatial filter — represented 
by the fiber — as a function of telescope diameter. This 
combines essentially the work of Shellan [34J — which de- 
scribes the on-axis intensity after AO correction — with 
the specific spatial filtering of a fiber, which is roughly of 
Gaussian shape and therefore de-emphasizes the role of 
higher modes in Zernike expansions. 

II. OPTICAL SETUP AND MODEL 
A. Fibre Set 

Step-index fibers are modeled with core radius a, re- 
fractive indices n c in the core and n in the cladding, and 
the same numerical aperture 

a m = 1.4 = V^c + n2 (!) 

for all bands. The difference in n c — n is kept at 
0.36 percent — copied from a specification of the Corn- 
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ing SMF-28e Photonic Fiber. This is rather hypothetical 
since we do not look at the absorption characteristics of 
materials in the infrared [H, [l8| . The cutoff wavelengths 
A c considered here are derived from the normalized fre- 
quency v w 2.405 [HI 

v = aka m , (2) 

where k = 2n/X is the momentum number, which leads 
to the fiber specifications of Table U 



TABLE I: Choice of fiber geometries, cutoff wavelengths A c 
and core radii a in infrared bands. 



band 


A c (/xm) 


a (/an) 


n c 


n 


J 


1.13 


3.089 


1.65437 


1.64843 


H 


1.36 


3.718 


1.65437 


1.64843 


K 


1.87 


5.112 


1.65437 


1.64843 



The two-dimensional distribution of the electric field 
at the fiber's entrance is computed by matching the two 
Bessel Functions in the core and in the cladding at each 
individual wavelength A [III HfiJ • The electric field func- 
tions f(p) depend only on the distance p to the fiber 
axis (and parametrically on A, a, n c and n), since no 
azimuthal dependence is left in the monomode regime. 
A numerical Hankel transform 0, @, d, d E H& H| 
transforms this into the far field F(r) 

/>OC 

F(r)oc / dppf(p)Upka) (3) 
Jo 

in the pupil plane. (No Gaussian approximation to the 
far field [3l|, [44j is introduced.) The integration over the 
azimuthal angle over the fiber's cut has already been per- 
formed in the Fraunhofer approximation, and has been 
condensed to the Bessel Function Jo. r is the radial co- 
ordinate in the pupil plane, a < a m the angle from the 
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fiber axis to this point in the pupil plane. The factor p 
is the Jacobian from the introduction of circular coordi- 
nates in the plane of the fiber's front face. The far field 
does not depend on the azimuthal angle 9 in the pupil 
plane. 



B. Phase Screens 

The electric field in the exit pupil of the telescope is 
written as a two-dimensional phase screen over the radial 
coordinate r and azimuthal coordinate 6 as 



E(r) 



(4) 



Amplitude variations — scintillation as opposed to phase 
variations, or imaging characteristics [H, [H, [U [3(| 
H3| — are not studied here, so the phases <p(r) are kept 
real-valued. Since we shall look only at coupling coeffi- 
cients, the modulus of E(r) is arbitrarily normalized to 
unity. 

The phase screens have been generated with a Kol- 
mogorov spectrum by synthesizing Karhunen-LoeyefKL) 
basis functions as described in the literature 0, [22|, 129 . 
S3: 



p.q 



We use normalized azimuthal basis functions 



e ? (0) 




cos q6 
sin qO 



(5) 



(6) 



(7) 



is Neumann's factor. If the radial basis functions are 
normalized according to 



where 




D/2 



rK^(r)K { /\r)dr = 5 pp ,, 



(8) 



the variance of the expansion coefficients is 

N^g p , q = D 2 [D/r»{\)fl*Bl ql (9) 

where B 2 are the eigenvalues of a reduced KL equation 
Q. Note that there is some arbitrariness in distributing 
factors here: the factor D 2 in @ might be absorbed in a 
renormalization of K p q \ and the factor ^/e 9 /(27r) of §6§ 

could also be dispersed over and/or B. 

The turbulent atmosphere is represented by a Kol- 
mogorov power-law of the phase structure function — we 
do not discuss the validity of this ansatz from any fun- 
damental or experimental point of view. The radial ba- 
sis functions are generated numerically by solving the 



symmetrized integral equations for the eigen-modes K p q ^ ; 
Zernike polynomials 0, [H, HH, HH, 13] or polynomials fits 
J4] have not been employed. 

The Fried parameters tq that are quoted here are those 
implicitly measured at 0.5 /J,m and were actually scaled 
with @ 



r oW — r 0|0.5/jm ( 



A 



0.5/^m 



6/5 



(10) 



to the infrared wavelength A to generate the kernel of 
the KL integral equations. The radial functions K p q ^ are 
re-generated for each instance of the ratio D/ro(A). 

An Adaptive Optics (AO) correction parameter (de- 
gree) c is introduced which assumes that some set of 
low-order basis functions with largest eigenvalues B p ^ q is 
discarded while building the full phase screen, c > 2 
means tip and tilt are removed from each individual 
phase screen, represented by the first line in [43|, Table 
III]; c > 5 means correction through refocus and astig- 
matism and removal of the first three lines in [43|, Table 
III], and c = 14 assumes correction of the modes of the 
first eight lines of [HI, Tabic III]. The calculations were 
done on a set of 75 basis functions, sufficiently large in 
comparison to these low-order corrections. 

The survival of speckles is demonstrated by coupling 
800 phase screens into a virtual photometric channel as 
if one would project a single telescope's input, directly 
with the coupling lens onto the fiber head. 



C. Coupling 

The intensity coupling efficiency is computed from nu- 
merical evaluation of overlap integrals over the circular 

pupil [HEI, 



.4 



\f NA E(r)F(r)d 2 



f NA E*(r)E(r)d 2 r f NA F*(r)d 2 r 



(11) 



The phase screen samples are created by generation of 
independent phase screens from uncorrelated Gaussian 
random numbers for the expansion coefficients g p , q ] in 
that respect no time scale is needed to define the transit 
time from one sample of the Kolmogorov statistics to 
another (l5l |40T ] ; there is no such parameter as the wind 
velocity of Taylor screens [l(| ■ 

Geometric imperfections like fiber misalignment (39l . 
l4lj ] or splicing are not incorporated, nor Fresnel reflec- 
tion losses at the fiber front end or a central obscuration 
(shadow of a secondary mirror) [l3|. 

For the Unit Telescopes of the Very Large Telescope 
Interferometer equipped with AO this has been studied 
in Fclkcl's thesis [3] • 

For weak turbulence, equation ([4]) can be expanded 
into its Taylor series 1 + iip — <p 2 /2 ■ ■ -. This decomposes 
the overlap integral in the numerator of pip into a con- 
stant (which represents the limit of optimum coupling 
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efficiency, ro — > oo), linear terms of azimuthal modes 
proportional to sin 6 or cos 9 which do not contribute be- 
cause the integral vanishes, plus quadratic terms which 
are quadratic in g Piq Bp q . The assessment of Ruilier and 
Fried that the energy coupling is represented by these 
terms remains basically correct, although the smearing 
with F(r) diminishes the influence of terms of larger ra- 
dial nodal numbers p. 



III. PHASE SCREEN STATISTICS 

Statistics over 800 phase screens have been compiled 
for sub-average (ro = 10 cm) and better-than-average 
(ro = 20 cm) seeing conditions [l(| HH for telescopes of 
D = 1 m (Figures [THB]) and 1.8 m in diameter (Figures 
[7 HT2]) at AO corrections levels of c — 2, 5 and 14. 

In each band, illumination by three different wave- 
lengths A is studied, and the figures that follow are split 
into three panels. Standard diffraction theory shows how 
the blue components of spectra are enhanced for cou- 
pling into a waveguide of fixed geometry [27j], but there 
is a counter-effect through the smoothing of the phase 
structure functions with (fTU|) . Here, the net effect is a 
higher coupling efficiency for the red end of the bands, 
in accordance with earlier results [H, Fig. 1]. This is 
not necessarily the full truth since the phase screen basis 
functions are calculated for each individual wave- 

length as a function of its phase structure function; as 
our model leaves a number c of these aside, we are im- 
plicitly presuming that the AO performs on the same 
level at all wavelengths within a band, and this might 
not be realistic. 

Tables [TT] and IIIII summarize the cumulative distribu- 
tion function of each statistics by a triplet of values in- 
dicating the median (50 % percentile) and the offsets 
from there to the 84.1 % and 15.9 % percentiles, equiv- 
alent to providing error bars on a la level. The nota- 
tion is A 50 + A 8 4.i — A15.9 where 50 percent of the cou- 
pling efficiencies are smaller than A50, 15.9 percent are 
smaller than A50 — A15.9, and 15.9 percent are larger than 
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FIG. 1: J-band. D = 1 m. r (500ran) = 0.1 m. 




FIG. 2: J-band. D = 1 m. r (500nm) = 0.2 m. 



4 





FIG. 3: H-band. D = 1 m. r (500nm) = 0.1 m. 



FIG. 5: K-band. D = 1 m. r (500nm) = 0.1 m. 
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FIG. 4: H-band. D = 1 m. r (500nro) = 0.2 m. 



FIG. 6: K-band. D = 1 m. r (500nm) = 0.2 m. 
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TABLE II: Summary of percentiles of coupling efficiencies 
for a telescope diameter of D = 1 m. 



band 


i o ^uuu 11111 j 
l m J 


Fio- 
r ig. 


(/illl) 




/l +-4 84 .i 
A50 -A 15 .9 


T 
J 


n t 
U. 1 


m 


1.1 i 


Z 


n or + 0.13 

u.oo_ 14 


7 

J 


U. 1 


m 


1 17 
1.1/ 





n C4+0.07 
U-O4_ 10 


J 


U. 1 


"i i 


1 OK 
l.ZO 


z 


r, OQ+013 

u -° 3 -0.14 


■ J 


U. 1 


m 


1 9^ 
l.ZO 





n c- 7 +0.06 
u - J '-0.09 


T 

■J 


o i 

U. 1 


i 


1 ^9 
X .OZ 


9 
z 


n 49+0.12 


J 


o i 

U. 1 


m 


1 *X9 
X .OZ 





n cq+0.06 

u.t>y_ 09 


T 

J 


n 9 
u.z 


9 


11V 
1.1 ( 


9 

z 




J 


u.z 


9 


11V 
1.1 ( 





n 71 +0.03 
U -' i -0.04 


T 
J 


o 9 
u.z 




1 9^ 
l.ZO 


9 
z 




T 
■) 


u.z 


- 


1 9^ 
l.ZO 





n 7o +0.02 
U-^2_ 04 


T 

J 


n 9 
u.z 




1 ^9 
1 .OZ 


9 
z 


n cc+0.05 

U.DO_ 08 


■ J 


u.z 


9 


1 Q9 
l.oZ 





n yo+0.02 
u - ' °-0.03 


TJ 
xx 


n i 

U. X 


Q 
■5 


1 4S 

X .4:0 


9 


n 48+ - 11 


tj 
XI 


o i 

U. 1 


•) 
■> 


1 /IS 





n cq+0.05 


TJ 

xi 


n i 

U. 1 


•V 
■J 


X.DO 


9 
z 


n rq+0.09 


TJ 


U. 1 


121 


l.DO 







TJ 
Xl 


o i 

U. 1 


•3 


1 7(\ 
1. 1 O 


9 
z 


n cc:+0.08 
u - dJ -0.12 


XI 


U. 1 


•3 
'> 


1 Vfi 
1.(0 







XI 


n 9 
u.z 


4 


1 /IS 
1 .^Xo 


9 
z 




TJ 
X X 


n 9 
u.z 


4 


1 /IS 
1 .^xo 





n 74+0.02 
U -' 4 -0.03 


TJ 
XI 


n 9 
u.z 




l.DO 


9 
z 


n 7n +0.04 


TJ 
XX 


o 9 
u.z 


4 


l.DO 





n 7 c:+o.oi 


TJ 
XX 


n 9 
u.z 


A 
4 


1 7f\ 
1. 1 O 


9 
Z 


n 71 +0.03 


TJ 
XI 


n 9 
u.z 




1 7f\ 





n 7 c+o.oi 

u - ' J -0.02 


T^ 
XV 


1 

U. 1 


r. 

a 


9 09 
z.uz 


9 
z 


n ri+ 07 


T^ 
XV 


1 

U. 1 


e 
a 


9 09 
Z.UZ 





n 7n +0.03 


K 
XV 


n i 

U. 1 


r. 

2 


9 9^ 
z.zo 


9 
z 


n C4+O.O6 


XV 


n i 

U. 1 


r. 

12J 


9 9^ 
Z.ZO 





n 7 9+0.02 
U -' Z -0.04 


T^ 
XV 


o i 

U. 1 


r. 

a 


9/11 
Z.^l 


9 
z 


n fic+0.05 


XV 


n i 

U. X 


a 


2.41 


n. 
U 


n 7 9+0.02 
u - ' ^-0.03 


K 


0.2 


6 


2.02 


2 


n 70+0. 03 
u - ' °-0.04 


K 


0.2 


6 


2.02 


5 


u -' '-0.02 


K 


0.2 


6 


2.25 


2 


n 74+0-02 

u - ' 4 -0.03 


K 


0.2 


m 


2.25 


5 


77+0- 01 
u - ' ' -0.01 


K 


0.2 


6 


2.41 


2 


n 7^+ 02 

u - ' J -0.03 


K 


0.2 


6 


2.41 


5 


77+0- 01 
u -' '-0.01 




FIG. 7: J-band. D = 1.8 m. r (500nm) = 0.1 m. 




FIG. 8: J-band. D = 1.8 m. r (500nm) = 0.2 m. 
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FIG. 9: H-band. D = 1.8 m. r (500nm) = 0.1 m. 



FIG. 11: K-band. D = 1.8 m. r (500nm) = 0.1 m. 
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FIG. 10: H-band. D = 1.8 m. r (500nm) = 0.2 m. 



FIG. 12: K-band. D = 1.8 m. r (500nm) = 0.2 m. 



7 



TABLE III: 
D = 1.8 m. 



Summary of percentiles of coupling efficiencies, 
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IV. SUMMARY 

A compact view on the median coupling efficiencies of 
these two tables is given in Figure [T3l The values on the 
curve with the red crosses, c = 2, are slightly more opti- 
mistic than those of Shaklan-Roddier [33|, Fig. 2] which 
one may attribute to underestimation of efficiencies by 
the Gaussian approximation. 

By dividing the median through the noise introduced 
by the Kolmogorov fluctuation in the phases, we ob- 




c=2 + 
c=5 
C=14 * 



++, ^x**. 



2 3 4 5 

D/r„ 




FIG. 14: The signal-to-noise ratio represents the median cou- 
pling efficiency divided by the statistical width of the distri- 
butions (la error bars). Both plots represent the same data 
set, one marker for each of the curves in Figures [T1 412I the 
upper on linear and the lower on logarithmic axis scales. 
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tain A 50 /(A S 4,i — Ai 5 , 9 ) as a signal-to-noise ratio (SNR) 
for each combination of fiber geometry, wavelength, tele- 
scope diameter, Fried parameter and AO correction. Fig- 
ure H3] shows these twice. The aim of the doubly loga- 
rithmic representation is to demonstrate that the SNR 
can be well fitted by a power law cx (D/ro)~ 5 ^ 3 with a 
prefactor depending only on the AO correction level c 
|32| | . This is not unexpected because this power has been 
an input to equation Q, and for a pinhole type of spa- 



tial filter this remains essentially unharmed as reasoned 
in Sect. mO[li|. 

Reduced coupling efficiencies as a function of the num- 
ber of speckles and degree of AO correction have been 
well reported by Shaklan and Roddier. We have illus- 
trated that in addition this reduction in photometric sig- 
nals is accompanied by wider variances of the expected 
coupling efficiency, which leads to the need of longer in- 
tegration times during observations. 
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